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Abstract 

A discrete model is introduced to account for the time-periodic oscillations 
of the photocurrent in a superlattice observed by Kwok et al, in an undoped 
40 period AlAs/GaAs superlattice. Basic ingredients are an effective negative 
differential resistance due to the sequential resonant tunneling of the pho- 
toexcited carriers through the potential barriers, and a rate equation for the 
holes that incorporates photogeneration and recombination. The photoex- 
citing laser acts as a damping factor ending the oscillations when its power 
is large enough. The model explains: (i) the known oscillatory static TV 
characteristic curve through the formation of a domain wall connecting high 
and low electric field domains, and (ii) the photocurrent and photolumines- 
cence time-dependent oscillations after the domain wall is formed. In our 
model, they arise from the combined motion of the wall and the shift of the 
values of the electric field at the domains. Up to a certain value of the pho- 
toexcitation, the non-uniform field profile with two domains turns out to be 
metastable: after the photocurrent oscillations have ceased, the field profile 
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slowly relaxes toward the uniform stationary solution (which is reached on a 
much longer time scale) . Multiple stability of stationary states and hysteresis 
are also found. An interpretation of the oscillations in the photoluminescence 
spectrum is also given. 

Key words: current instabilities; electric field domains; metastability; dynam- 
ical effects in superlattices. 
Pacs numbers: 72.20.Ht., 05.45+b, 85.30.Fg 
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I. INTRODUCTION 



In recent experiments by Kwok et al. the time-dependent transport properties of pho- 
toexcited undoped superlattices were investigated |^. In a typical example, an undoped 40 
period AlAs/GaAs superlattice was mounted in a p-i-n diode and continuously illuminated 
by laser light at 4K. When the laser power was in a certain interval and the applied dc volt- 
age bias was large enough, damped time-dependent oscillations of the photocurrent (PC) 
and the peaks in the photoluminescence (PL) spectrum were observed. The applied fields 
in the experiment [|T| are high, the potential barriers are wide (between 30 and 40 A), and 
the minibands correspondingly narrow so that the coherence length is comparable or smaller 
than the width of one quantum well. Then the quantum wells in the superlattice are weakly 
coupled and formation of electric field domains may appear [^^. The subject of formation 
and expansion of electric field domains in (doped) superlattices is an old one and it goes back 
to the conductance measurements of Esaki and Chang in 1974. Typically the dominant 
transport process in each domain is resonant tunneling between different subbands of the 
wells belonging to the domain. 

Time-dependent oscillations of the PC in undoped illuminated superlattices have also 
been observed by Le-Person et al. ^ for quite different superlattices where miniband trans- 
port is prevalent (narrow barriers, so that the quantum wells are not weakly coupled. In 
these conditions it is not clear whether formation and propagation of electric field domains 
are relevant). These measurements were interpreted in terms of Gunn dipole domains of a 
classical drift-diffusion model. 

In experiments with weakly coupled superlattices at high electric fields, the dominant 
transport mechanism is resonant tunneling: It follows from Kazarinov and Suris's tight bind- 
ing calculation that when the Wannier-Stark level splitting is smaller than the distance 
between the two lowest subbands of each quantum well of the superlattice, the electron is 
localized and the current through the structure is small. When both quantities are com- 
parable, there is a peak in the current-field curve due to resonant tunneling between the 
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first subband, ei, of one well and the second subband, 62, of the next one (let us denote 
by Ci —>■ 62 this kind of resonant tunneling). Further peaks correspond to other resonant 
transitions between different energy subbands of the wells. As a result of resonant tunnel- 
ing, electrons populate higher subbands of the wells and PL measurements may be used to 
probe the change of electron population in higher subbands. [Q,^ In superlattices with wide 
wells the resonant tunneling process is sequential, which means that after a tunneling event 
the electron decay very fast to a lower (e.g., the first) subband due to scattering. For 



the superlattices described above, the intersubband scattering time is of the order of 1 ps 
which is very small compared to the typical tunneling time of 100 ps. When electric field 
domains are formed, the wells belonging to each domain undergo similar sequential resonant 
tunneling (SRT) events. In Ref. and previously in Ref. 0], three domains are observed 
for different bias and illumination: domain I where SRT is ei ei (the subbands of all 
the wells in this domain are aligned at zero electric field), domain II with ei 62 SRT, 
and domain III with 62 — > 63 or ei — > 63 SRT according to the sample. We have depicted 
schematically the SRT processes at a sample where domains II and III coexist in Fig. 
Time-dependent PC oscillations are observed as a transient toward a stationary state where 
domains II and III coexist. |^ 

In this paper we describe and analyze a simple discrete drift model that explains several 
features of the experiments reported in Ref. [|l|, particularly the stationary current- voltage 
characteristic curve and the time-dependent oscillations of the PC in the region of the 
parameter space where domains II and III coexist. Previous theoretical works on electric 
field domains in superlattices include Grahn et al's characterization of stationary states. 



Laikhtman's equivalent circuit model [|T2| and Korotkov et al's simulations based upon the 
analogy of a slim superlattice with a one-dimensional array of metallic tunnel junctions. 
13| We shall comment on these works in the Discussion. To motivate our model let us 



consider several important time scales for the superlattices object of our study. First of 
all, the characteristic time scale of the PC oscillations in Ref. [ffll is of 100 ns. The time 
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scale for carrier thermalization is of 0.1 ps while the carriers reach thermal equilibrium with 
the lattice after a time that ranges from 1 to 100 ps. [|T^ This means that in time scales 
of the order of nanoseconds (the experimental time scale), we may consider the holes and 
electrons to be at local equilibrium within each quantum well j at the lattice temperature, 
and with given values of their densities, pj and fij. The process of reaching a stationary 
state might be seen as the attempt of reaching a "global equilibrium" starting from "local 
equilibrium" through tunneling processes that communicate different quantum wells, self- 
consistency of the electric field and scattering and interband processes. In this spirit we 
consider the quantum wells as entities characterized by average values of the electric field, 
Ej for the j-th well, and of the densities of holes and electrons, pj and nj, respectively. We 
then propose the following transport equations to describe the dynamics of the superlattice: 

el ^ 

Ej - Ej_i = — {uj -pj), (1) 

eS + ei)(E,)n, = J, (2) 
dpi ^ ^ ^ ^ 

—^ = ^-rnjPj, (3) 

where j = 1, . . . , N. In this model: 

1. Equation (|l]) is the discrete Poisson equation relating the field at two adjacent wells 
and the electric charge; e is the electron charge, / is the length of one period of the 
superlattice, and e is an average permittivity. Equation (|l|) may be considered as the 
average of the usual Poisson equation over one period of the superlattice. Since the 
heterostructure corresponds to an intrinsic semiconductor, we do not include charged 
impurities in (|l|). Poisson's equation for the j = 1 well contains the field at the 
"zeroth" well, which corresponds to the doped semiconductor before the superlattice 
and it will be dealt with below. 

2. Equation (|^) is Ampere's law at the j-th well establishing that the total current density 
J is the sum of the displacement current and the electron flux. We consider that the 
holes are heavy and their contribution to the current is negligible. [0 The electron 
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flux JjL,j-^-i = ^v{Ej)nj is due to sequential resonant tunneling between subbands of 
neighboring quantum wells as said above. By writing this expression we are assuming 
that the probability of tunneling from the well j to the well j ' + 1 is proportional to the 
number of electrons in the well j, and we are ignoring the small "reverse tunneling" 
from the well j ' + 1 to the well j. The probability of tunneling is larger for fields 
corresponding to SRT (for which the subband ei of j is aligned with 62 of j + 1, or 62 
of j with 63 of j + 1, etc). This implies that the function v{E) (with the dimensions 
of a velocity) has peaks at the corresponding values of the electric field. A crucial 
assumption of our model is that the electron "velocity" corresponds to the static 
current- volt age characteristic of the superlattice at low laser power (small 7; see (^ 
below). We take this curve as a datum of our model and argue that our results do not 
depend qualitatively of the exact shape of v{E). Since we are interested in the analysis 
of the region where domains II (ei 62 SRT) and domains III (e2 63 or ei 63 
SRT) coexist, we have used a velocity curve with two local maxima centered in the 
same field values as seen in the experiments. The result is Fig. ^. To explore the 
region where domains I (ei —>■ ei SRT) and II coexist, it is convenient to add another 
(smaller) peak to the velocity curve at zero electric field. The modifications will be 
reported elsewhere. 

3. Equation (^) is the rate equation for the hole concentration. We consider only processes 
of photogeneration of electron-hole pairs and recombination. The photogeneration 
term, 7, is proportional to the power of the laser (Ref. |16|, chapter 12; see also 
Delalande's paper in Ref. jl^), which we assume (for simplicity) that it illuminates 
uniformly the superlattice. The recombination coefficient is assumed to be constant. 

Notice that the rate equation for the electrons containing the photogeneration, recombi- 
nation and SRT processes is: 

^=7-rnjPj + l{v{Ej_i) nj_i - v{Ej) fij}. (4) 
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Here the SRT current coming from well j — 1 to well j brings electrons to the j-th well at a 
rate v{Ej_i) hj_i/l, while the SRT current going from well j to well j + 1 detracts electrons 
from the j-th well. This equation is equivalent to (Q), as it can be shown by differentiating 
(I) and using and (|). 

In the 3A^ Equations there are 3A^+2 unknowns, J, Eq, Ej, hj,pj, with j = 1, . . . ,N. 
One additional equation is the bias condition: 

1 A ~ $ 

T7E^i = — - (5) 

Here $ is the difference between the applied voltage and the built-in potential due to the 
doped regions outside the superlattice (1.5 Volts in Ref. [|I|). is the average applied electric 
field on the superlattice, which we will henceforth call the bias. The missing condition is a 
boundary condition for the field at the zeroth well, Eq, {before the superlattice). We do not 
have direct experimental evidence for what Eo should be. Thus our choice for Eq has to be 
validated a posteriori by comparing the results of our analysis with experiments and with 
the consequences of a different choice. We shall use throughout this paper the following 
boundary condition: 

Eoii) = E,{i). (6) 

This condition does not allow a charge build-up at the first well and it is compatible with 
the experimental observation that the field is the same for all the wells at low laser power: 

for all j. This equation says that the effective electron velocity is (except for constant scale 
factors) the same as the static I-V characteristic curve of the superlattice at low laser power, 
which allows us inferring the form of v{E) from experimental data. Thus our boundary 
condition (^ forces a uniform field profile (the same Ej for all j) to be a solution of our 
model for any laser power. As we shall see later, this uniform solution is stable for low 
laser power whereas the model subject to the same boundary condition (P) may have stable 



non- uniform field profiles (with domains) for large enough laser power. Further discussion 
of the boundary condition is postponed to Section VI. For the velocity curve v{E), we 
have used a variety of phenomenological curves with two peaks and a negative differential 
resistance (NDR) region of negative slope between them. We could have also obtained the 
velocity curve from theories of static I-V curves for double barrier structures, but given the 
qualitative insensitivity of our results to variations of the velocity curve (as we will explain 
later), we have stuck to curves with easy analytical expressions for the calculations presented 
here. 

The boundary condition (^ allows for electric field domains Ej, j = 1, . . . ,N made of 
two constant values of the electric field and a domain wall connecting them, as we shall 
discuss below. The fact that our superlattice is embedded in a p-i-n diode implies that 
domains with low field at x = and high field at x = L are preferred to domains with 
the symmetric configuration (see Eq. ([l|)). Thus we shall only consider initial conditions 
that favor this kind of electric field domain and we will not mention the domains that start 
at X = with high electric field, even though they may be possible with our boundary 
condition. When the superlattice forms part of a different diode, e.g. n-i-n, the two kind of 
domains are allowed and richer structures may arise according to the initial condition used 
in the simulations. 

It is convenient to render dimensionless the equations of our model before we proceed to 
their analysis. For this, we adopt as the unit of electric field that at the first maximum of 
the velocity curve, Em — 10^ V/cm. Then (|1|) yields the typical unit of electron density 

fi = !:^^ioi7cm-3, (8) 

el 

for / ^ 100 A. There are two important time scales in our model. From (|^) we see that the 
electrons employ a time = 1/v{Em) — 0.4 ns to tunnel from one well to the next one of 
the superlattice. From (|^), we find the time scale in which the hole density varies 

— = fh. (9) 

Tp 
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Thus Tp is essentially given by the recombination coefficient r. f is a decreasing function of 
the electric field because the overlap between the electron and hole wave functions decreases 
as the field increases. This accounts for a variation of Tp from subnanoseconds to tens 



of ns. ||T8| Here we adopt Tp = 10 ns which is reasonable in view of the experimental values 
reported in Ref. We measure the time in units of Tp which together with the previously 
mentioned units complete the following definition of dimensionless variables: 
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M 



Pi 



n 



~ 5 
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t 



t = — = rnt, J 
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Inserting (ITOlO) into we find the dimensionless system 



(10) 



(11) 



dE- 

(3^ + v{E,)n, = J, 



dpj 
~dt 

1 ^ 
- TEj 



Eq = El. 



(12) 
(13) 
(14) 

(15) 
(16) 



Here /3 goes from 0.01 to 1 and and 7 are the dimensionless control parameters. These 
equations are to be solved with initial conditions for the fields, -Ej(O), and the hole concen- 
trations, Pj(0), compatible with the bias (|T5|) and the boundary condition (plGf). The initial 
conditions for the electron density, ^^(0), then follow from (|12|). 



II. STEADY STATES 



For given values of 7 and J there are stationary solutions of the equations ([T^-|T^ 
compatible with the boundary condition (|16]). Using equation (|l^), we find the corresponding 
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value of 0. From this we can obtain the curves J = J{<f)) appearing in Fig. ^ We have 
studied two different types of stationary solutions, namely uniform solutions and two-domain 
solutions. 



A. Uniform solutions 



By inserting Ej = E (constant) into (|T2|-^), we find 



E = <p (17) 
ri,=Pj = Vl U = l,...,N) (18) 
J=^v{<P) (19) 

Equation ([T9|) yield a static current- volt age characteristic curve which is the same as v{E) 
except for a scale factor. See Fig. |^. 

B. Non-uniform two-domain solutions 

The stationary versions of equations (|l^) and (|l^) can be solved for nj and Pj in terms of 
the electric field Ej and then inserted in (0). The result is the following discrete mapping 

E,_i = /(E,;7,J) (20) 

with 

/(E;7,J) ^E--^^+lv{E). (21) 

The mapping (pOD allows us to get all the electric field profiles {Ej, j = 1,...,A^} that 
satisfy the boundary condition Eq = Ei. From these profiles and the bias condition (|15|), 
we obtain the current- volt age characteristic curve of Fig. |^. Motivated by experiments, we 
are interested in stationary solutions with two electric field domains, i.e. Ej = El {El is 
the constant low field value) for j = 1, . . . , k; Ej = Eh {Eh is the constant high field value) 
for j = N — k + 1, . . . , N, whereas for the remaining quantum wells {j = k + 1, k) Ej is 
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an increasing function of j that takes values on the interval {El, Eh)- These wells constitute 
the domain wall that separates domains II and III. 

For these solutions to exist, a few requirements have to be fulfilled. First of all, the 
boundary condition ( [T6| ) together with (^) imply that Eq should be a fixed point of the 
discrete mapping (^), i.e., a solution of 

fiE;^,J) = E (22) 

(see Fig. §). Furthermore in order to have a non- uniform solution with domains, (^) should 



have more than one fixed point. After simplification, Eq. (^) becomes (|^) and therefore 
this requirement is satisfied for certain values of J described below. Let E^ be the fixed 
point of the mapping that we choose for Eq. With the v{E) profile depicted in Fig. |^, 
(^) may have one or three solutions depending on the values of J and 7. Let Em, Em' 
and Em be the positive electric field at the first and second positive maximum and at the 
first positive minimum of v{E), respectively (see Fig. ^). For < J < ^v{Em) and 
for ^v{Em) < J < ^J^v{Em'), ( P^ has only one solution on the interval {Q,Em')- For 
y/^v{Em) < J < ^v{Em), there are three solutions of ( p2D on the same interval. Thus 
we can find a fixed point E^ on the interval {0,Em') for < J < ^v{Em')- Secondly, a 
non-uniform profile should have k < N wells with fields equal to El and Ek+i > El- This 
means that the equation 

f{E;j,J)=EL, (23) 

should have at least a solution different from the fixed point, otherwise only the uniform 
solution described above is possible. This second condition holds only for 7 > 71 ^ 5.2 x 10~^ 
(see Section |11I[ ). If one of the two lowest fixed points of (^) satisfies these two conditions, 
we can select it as the value of the low field domain El and let the k first wells to have this 
value for their electric fields {Ej = El, Vj = 1, . . . ,k). Then one may choose as Ek+i another 
solution of ( P5D different from El- The p-i-n character of the diode prescribes that this is 
an acceptable solution provided the electric field Ek+i is higher than El- Then (^) implies 
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that Ej increases with j and approaches another fixed point, which is the value Eh of the 
high field domain. As Fig. ^ shows, for each fixed 7 > 7i both conditions are satisfied only 
for values of J on a certain interval and for some of the fixed points of the mapping. The 
current- applied bias {J — (j)) characteristic curves obtained with the previously described 
method are depicted in Fig. ^ for 7 = 0.016. As explained before, the continuous curve 
(together with the short-dashed part in the NDR region) is ^/7f (0), with v{E) depicted 
in Fig. 1^. It corresponds to the uniform solution, obtained when all the Ej stay on one of 
the fixed points (^). The J-cf) curves for the non-uniform stationary solutions with two 
domains explained above appear for (j) on an interval that ranges from a value slightly below 
Em to a value slightly below Em' in Fig. ^. They exist for 7 > 71 and J on a subinterval of 
(V^'^(-^"i)' V^'^i^M)) whose width increases with 7. 

Let us now describe several important non-uniform stationary profiles of the electric field. 
Consider in the first place the curves formed by small triangles in Fig. 0. They correspond 
to non-uniform solutions with Ej = E^ (the first fixed point of the discrete mapping on the 
first branch of v{E)) for I < j < k (the number k is shown on top of several such curves 
in Fig. 1^). For j > k, Ej "jumps" to the third branch of f{E;'y, J) and it approaches the 
third fixed point of the discrete mapping (1 — > 3 jumps). See these iterations of the discrete 
mapping in Fig. §(c) for a typical value of J. As increases and k decreases, each branch in 
Fig. ^ corresponds to the same steady field profile but with the domain wall displaced one 
step backward (towards x = 0), which means that a new quantum well has lost the low field 
value. In the main body of Fig. ^, only one out of each three branches of stationary profiles of 
this type are depicted. The inset shows all these solution branches for k between 17 and 25. 
Since they turn out to be linearly stable (see the Appendix), multistability between different 
stationary profiles corresponding to the same voltage (p is possible. As 7 grows, the width of 
the J interval for which these solutions exist increases, and so does the length of the triangle 
curves in Fig. |^. Then coexistence of more than two stable stationary filed profiles for the 
same is possible. This gives rise to larger hysteresis cycles and memory effects. When do 
the stationary field profiles with 1 — 3 jumps exist? As the sequence of Figs. §(a) to §(c) 
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shows, the values of 7 and J should allow the local minimum of f{E] 7, J) to descend below 
the first fixed point of the mapping. This is possible for 7 > 71 and J G {Ji{^), y/jv{EM)), 
where the minimum of f{E] 7, J) coincides with the first fixed point at the current ^1(7). 
The branches of stationary solutions with 1^3 jumps end at J = ^v{Em) when the fixed 
points on the first and second branches of the discrete mapping coalesce and disappear. 

Consider now the curves with small squares (and their long-dashed continuations) in Fig. 
^. They correspond to stationary solutions with Ej = El, which is now the second fixed 
point of the discrete mapping, for the first k quantum wells. For j > k, Ej jumps to the 
third branch of f{E; 7, J) (2 — > 3 jumps). The short-dashed lines in Fig. ^(c) represent the 
iterations of the discrete mapping leading to one of these solutions. For a given 7, for which 
range of J do the stationary field profiles with 2^3 jumps exist? We shall argue (and have 
checked numerically) that the profiles with 2^3 jumps exist for J G (^2(7), y/l '^{Em)), 
where J2(7) is determined from (|T9|) and (^4]) below. From Fig.^b) we see that these profiles 
appear at a value of J = ^2(7) where the slope of f{E;j, J) at the second fixed point of 
the mapping is zero. We observe in Fig. ^ that all 2 ^ 3 solution branches taper down to 
a point 02(7) for J = J2{l)- We find 02 (7) and J = ^2(7) by inserting (|IU|) into the 



condition that the slope of f{E; 7, J) be zero: 

v + 2^v' = 0. (24) 

This equation has one solution corresponding to the minimum of f{E; 7, J) which determines 
02(7)- Substituting 02(7) into ([T9|) we obtain ^2(7)- The smallest possible 7 for which (^4]) 
has solutions is 7 = 71 ~ 5.2 x 10^^: Non-uniform stationary profiles exist for 7 > 71. As 
we indicate in the Appendix, the curve ( ^4]) bounds the region inside which the uniform 
stationary solution is linearly unstable. As in the case of the stationary field profiles with 
1^3 jumps, only one out of three branches corresponding to profiles with 2^3 jumps 
is depicted in Fig.|^. Clearly many more stationary solutions can be constructed with the 
help of the discrete mapping, and we have depicted only the more significant ones in Fig. |^. 
We want to point out that all the stationary branches are connected: they do not disappear 
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as it could be wrongly inferred from our Fig. 0. For example, in Fig. H(c) we have shown 
(long-dashed line) the mapping corresponding to a stationary field profile different from the 
pure 1^3 (continuous lines) and 2^3 (short-dashed lines) jumps. For these (unstable) 
solutions, -Efc+i jumps to the second branch of f{E; 7, J) (instead of jumping directly to the 
third branch) before continuing toward the third fixed point of the mapping. Clearly the 
branches corresponding to these solutions coalesce with the 1^3 branches at J = </i(7) 
and explain the disappearance thereof. 

Remark. We have found stable non-uniform stationary states having two domains with 
electric fields El and Eh separated by a domain wall. These fields are: 

• El < Em (curves formed by small triangles in Fig.||) or El G {Em, Em) (curves formed 
by small squares in Fig.|^). 

• Eh G {Em-, Em')- 

Thus while the first domain may have an electric field below or above the SRT value Em-, 
the field at the second domain is always below the second SRT value Em'- Notice that 
this happens independently of the applied bias for G {Em,Em')- That the field at the 
domains does not necessarily coincide with the field for which there is resonant tunneling 



is in agreement with experimental observations, and previous theories, |jT2[. For applied 
bias on {Em, Em'), the average of the stationary current at the branches representing non- 
uniform states is independent of the bias, which explains the "fiat" shape in Fig. |^. The 
experimental data present both fiat plateaus and rising profiles in the I-V characteristic 
curves of the samples. 



III. LASER POWER VS. APPLIED BIAS "PHASE DIAGRAM" 

With the aim of interpreting the experimental results, we have to delimit the intervals of 
the parameters 7 and (p where stable solutions of the types described in the preceding Section 
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may be found. A discussion of their linear stability is postponed to the Appendix, but this 
information together with that provided by the discrete mapping of the the preceding Section 
allows us to draw the 7 — "phase diagram" of Fig. ^. The continuous line in Fig. ^ bounds 
the region in which stable non-uniform solutions with domains exist. The flat bottom of 
this region is at 7 = 71, the photoexcitation below which non-uniform stationary states do 
not exist. The continuous line in Fig. |^ has been obtained numerically by the following 
procedure: 

From Fig. ^ we see that the branches of stable non-uniform solutions that yield the 
minimum and maximum values of (j) are the two extreme small triangles. We have used the 
discrete mapping (EH) to find with precision those extreme values of as functions of 7. For 



a given 7, the range of J for which the discrete mapping (pO]) allows for stable non-uniform 
solutions of any type (with 1 — 3 or 2 ^ 3 jumps) has been determined following the 
procedure explained in Section II. Notice that the minimum value of (in the continuous 
line of the phase diagram Fig. |^) is determined by taking the smallest J = ^1(7) for which 
it is possible to find a solution of ( pil| ) and ([TB|) having 39 wells with field at the lowest 
fixed point and 1 well with a higher electric field. On the other side, the maximum value of 
on the continuous line of Fig. |^ is reached at the end of the last small triangle branch of 
Fig. ||. This value of is found by letting Eq = Ei = El be the second fixed point of (^Of) 
at J = ^v{Em) (the value of J for which the two first fixed points of (^) coincide) and 
then using ( ^0]) to find the other values Ej > El with j > 1. 

Inside the continuous line in Fig. ^ we have depicted with a dashed line the curve (p4|). 
As we explain in the Appendix, the uniform solution is linearly unstable inside this curve. 
Between both curves in Fig. ^, both the uniform solution and at least one branch (often 
several ones) of non-uniform solutions are linearly stable for a given value of 0. As discussed 
in Section II, we can easily visualize mult ist ability with the help of Fig. ^ 



15 



IV. SIMULATION 



From the experiments we have information about the steady state of the problem through 
the current-voltage characteristics, and about the dynamical aspects thereof through the 
time evolution of the PC and of the PL. In our model we have direct access to the time 
evolution of the PC, proportional to J{t), while we would need to supplement our model to 
be able to obtain the time evolution of the PL. We shall consider the evolution of the J{t) 
and of the electric field and charge profiles in this Section, while we will give qualitative 
arguments about the evolution of the PL in the next Section. 

In the previous stationary study we have found a very rich 7 — J phase diagram, showing 
multistability and hysteresis. We have seen that our model accounts in a very satisfactory 
way for the static properties observed experimentally. The picture will be substantially 
completed by describing the time evolution of solutions of our model, which will explain 
most of the dynamical experimental results associated with the PC To extract time- 
dependent magnitudes from our model we have solved numerically the system of first order 
nonlinear equations, by means of a standard fourth order Runge-Kutta method. For a fixed 
value of voltage and laser power we turn on the laser at t = 0. The initial condition for 



the fields and the charge carriers are equations ([T^ - |T^, except that we increase the charge 
at the first quantum well to mimic the effect of the laser, and then let the system evolve 
according to Equations (|l^ - |16D . 

The analysis of the steady states in Section || reveals that at low laser power 7 < 71 
5.2 X 10^^ in our rescaled parameters), there is only one steady solution: the uniform 
one, ([IsHITD, which is linearly stable (to be more precise, the same stands for the whole 
region outside the full line of the j — (f) phase diagram shown in Fig. |^). This result confirms 
the usual assumption that the PC at low power laser is just a rescaled plot of the velocity of 
the electrons versus the electric field , and supports our choice of boundary conditions. 
Above this threshold there is always a voltage interval (inside the full line in Fig. ^ where 
multiple steady states exist. These states have step- like electric field profiles as described 
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above and some of them are linearly stable (see Fig. These solutions may even coexist 
with the uniform one, which is unstable only inside the dashed line in Fig. |^. 

Nevertheless, from a dynamical point of view, the simulations reveal an interesting be- 
havior in the laser power range < 7 < 72 (— 8 x 10^'^) for on a subinterval of {Em, Em') 
(for 7 = 0.016, 1 < < 1.4; see Fig. |^). The situation may be described as follows: 

A. < 7 < 71 

There is only one steady state: that corresponding to an uniform profile of the electric 
field. The linear stability analysis reveals that this steady state is stable. However, if the 
disturbance is large enough, the numerical simulation shows that, after a short time, a non- 
steady step-like electric field profile is formed. The domain wall of this profile oscillates 
back and forth in time about a fixed value of the position. As it oscillates, the domain wall 
changes its width, so that it spreads out over several wells, which makes easier the possibility 
of experimental detection. The electric field at the two domains in the profile also oscillates 
in time for a few periods before the oscillations stop (see Fig. |^ to visualize this process). 
Then there follows a very slow modification of the step-like profile which ends up in the 
uniform steady state. We call the step-like electric field profile that remains after the PC 
oscillations cease, the metastable state or the metastable domain formation. In Fig. |^, we 
have plotted the PC versus time for moderate and long times (inset), showing the oscillations 
and the final decay of the metastable state. Notice that the oscillations are appreciable only 
for 7 > 2 X 10^^. Thus for < 7 < 2 x 10~^ there is overdamped formation of a step-like 
field profile followed by a slow relaxation towards the uniform steady state. 

B. 71 < 7 < 72 

There are true non-uniform steady states some of which are stable. The numerical 
simulations show a pattern analogous to that described above without the slow drift towards 
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the uniform state: after the PC oscillations cease, a non-uniform steady state with two 
domains of the electric field is reached. 

C. 7 > 72 

The damping is so strong that there are no oscillations of the PC: the step-like electric 
field profile is reached in a monotone fashion. 

These results are visualized in Fig. | in which the evolution of the current with time 
is represented for increasing values of 7 in a range that covers the three situations just 
described. The region of the phase diagram (Fig. |^) where oscillations - about either stable 
or metastable non-uniform solutions - can be observed has been represented by a shaded 
rectangle: 2 x 10~^ < 7 < 8 x 10~^ and 1.0 < < 1.4, approximately. Outside this region, 
the damping about uniform or non-uniform (meta)stable steady states is so strong that no 
current oscillations may be observed. In Fig. ^ we have plotted the profiles of the electric 
field and of the electron and hole densities corresponding to a maximum and a minimum 
of the PC during one oscillation. Notice that the electronic charge oscillates from one side 
of the domain wall to the other during a period of the PC oscillation [|l|. These results 
provide a clear picture of domain formation and posterior time evolution accounting for the 
time-dependent oscillations of the PC. 

V. AN INTERPRETATION OF THE PL MEASUREMENTS 

The PC discussed in the previous section gives an integrated information about our sam- 
ple, and is by no means a clear indication of domain formation. The experimental evidence 
of domain formation is the evolution in time of the PL spectrum which measures the light 
emitted by the sample in the electron hole recombination processes. The time resolved PL 
is a direct way to measure the time evolution of the carrier distribution. A correct descrip- 
tion of the PL spectra would require explicit consideration of recombination light emission, 
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scattering and tunneling times in the SRT process and of experimental resolution, which 
is out of the scope of the present work. However, we would like to supplement our model 
with a few assumptions so as to get a qualitative interpretation of the PL measurements. 
To do this, we need to understand the processes involved on the basis of the experimental 
information [|^], which can be summarized as follows: 

1. For large enough laser power and applied voltage, two PL peaks appear when the laser 
is switched on, confirming the formation of electric field domains. Due to the Quantum 
Confined Stark Effect (QCSE) |^T| we can affirm that the peak that appears at lower 
energy corresponds to higher value of the electric field than that at higher energy. The 
intensity of each peak is proportional to the number of electron-hole recombinations. 

2. As the applied voltage increases the intensity of the high field peak (low energy) 
becomes greater but the position of the peaks remains approximately the same. 

3. The intensity of the peaks oscillates with time at the same frequency as the PC oscil- 
lations. 

4. The intensity of the low field peak (high energy) is greater for a maximum of the PC 
than for a minimum and the reverse occurs for the high field. 

The two first points indicate that the electric field value at each domain does not vary 
as the applied voltage increases: the domain wall just moves toward the p contact thereby 
increasing (decreasing) the size of the high (low) field domain. The third point shows that 
the electrons responsible for the PL are the same ones that carry the current across the 
sample. 

All these points are highlighted in the following very simple formula which relates the 
PL spectrum, PL{huj), to the carrier densities and the electric field at each well obtained 
from our model: 

N 

PL{huj) =fY.i^jPj^(^'^ - (1) 
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Here fij {pj) is the electron (hole) concentration at the well j, f is the recombination co- 
efficient assumed to be field independent, and Q{Ej) is the energy of the photon emitted 
in the recombination process at the j-th well. Due to the QCSE, this energy depends on 
Ej, and it can be easily calculated with a one- well model |22|. The result shows an almost 
linear behavior of Q{E) with slope 2elw (e is the electron charge and 1^ is the width of the 
well). In (|l|), A(x) is just a Lorentzian function that counts the contribution of ^{Ej) to the 
PL spectrum at energy hu with a finite resolution. In writing the expression ([l|) we have 
implicitly assumed that each well is at local equilibrium during the recombination process 
as explained in the Introduction. 

From Fig. ^ we see that a maximum of the PC occurs when the domain wall is almost 
vertical and the electric field at most of the wells takes on either the low {E^ — Em) or the 
high (Eh) field values, for which the tunneling probability is higher (cf. the curve v{E)). 
During one oscillation, the low field domain (which contributes to the high energy peak in 
the PL spectrum) is larger at the maximum of the PC than at its minimum. According to 
(|I|) and our previous discussion, this should imply an increase in the low field peak (high 
energy in the PL spectrum) for the maximum of the PC. When the PC takes on its minimum 
value within one oscillation, the domain wall is not so sharply defined and several quantum 
well have a value of the electric field for which the tunneling probability is small. 

The number of wells that are in the high field domain increases with the applied voltage, 
and therefore, the weight of the low energy peak (HF) increases with the voltage, refiecting 
the growth of the HF domain. In Fig. ^ we show the PL spectrum for a maximum and a 
minimum of the PC for the same values as Fig. ^ and ^. The result resembles the experi- 
mental spectrum but it is clear that in order to get quantitative agreement a more elaborate 
model of the sequential resonant tunneling processes and electron-hole recombination is 
needed. This is out of the scope of the present work. 
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VI. DISCUSSION 



We have analyzed a very simple model that addresses in a satisfactory way the transport 
aspects of the problem which are dominant at the time scales experimentally relevant (1 - 
100 ns |l|): sequential resonant tunneling between weakly coupled quantum wells, photogen- 
eration and recombination of electrons and holes, and self-consistency through the Poisson 
equation. An important aspect of our model is the discrete nature of its governing equations 
(Pl^rS). Recall that this crucial feature is imposed by the separation between the relevant 
time scales of the above said transport processes and those of intra- and intersubband scat- 
tering (cf. Section I). Our results exhibit a wide variety of phenomena in agreement with 
the experiments: multistability of the stationary solutions with its consequences of hys- 
teresis and memory effects @], time- dependent oscillations of the PC, domain and domain 
wall formation [|l|]. The metastability of the non- uniform stationary states predicted from 
our model is compatible with current experimental data, although longer observation times 
seem necessary for non-ambiguous confirmation |^. Notice that the non- uniqueness of the 
stationary solutions is a direct consequence of our equations and it can be easily explained 
with the help of the multivalued discrete mapping of Section |n|. The non-uniform stationary 
field profiles exist on the interval of applied bias {Em-, Em') which is significantly larger than 
the NDR region {Em, Em)- The region of multistability in the current- voltage characteristic 
curve is "flat" (the average current does not change with the applied bias, Fig.||). From our 
results it is clear that similar multistability due to coexistence between domains I (at zero 
field) and II (at field slightly below Em) @] may be obtained by adding a local maximum at 
zero field to the v{E) curve. A simple modification of our Poisson equation and the equation 
for the holes may allow the study of similar phenomena in doped superlattices, which will be 



considered elsewhere p3|. It is important to notice that our results hold for a generic curve 
v{E) of the form discussed in the paper and that we could fine-tune experimental data by 
changing the quantitative features of v{E) and the parameter values. For instance, we could 
increase the size of the region where PC oscillations appear (shaded rectangle in Fig. ^ 
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relative to {Em, Em') by displacing the local minimum Em towards Em'- As there would be 
PC oscillations about non-uniform steady states with larger biases, the domain wall during 
a period of these oscillations would be closer to x = than in Fig. |^. This would mean that 
we could obtain a larger domain III. In turn, this would cause a larger HF peak to appear 
in the PL spectrum of Fig. |TD| thereby improving the agreement with experimental results. 

Similarly the frequency of the PC time-dependent oscillations is a increasing function of 
/3: We could estimate the ratio between the tunneling and recombination times by fitting 
the numerically obtained values of the frequencies to the experiments. 

An important question, raised in the Introduction, is the dependance of our results with 
the boundary condition (^. Under this condition, the uniform stationary solution is stable 
at low laser power, which allowed us to identify the velocity curve with the I-V character- 
istic curve (except for unimportant constant scaling factors). At higher photoexcitation, 
the uniform stationary solution could become unstable in favor of non-uniform two-domain 
solutions which agrees with the usual experimental interpretation [Q. All these facts are of 
course indirect evidence in favor of our choice (^. Still one wonders whether other choices 
also yield results in agreement with the known facts. Without entering in detailed modeling 
of the contact and doped regions before the superlattice, one can consider other reasonable 
choices: it is plausible that modeling the region before the superlattice one obtains an in- 
creasing current- voltage curve (at least for not too large bias). Since Eq is the average field 
before the superlattice, one could derive a formula Eq = g{J) from that current- voltage 
curve and use it as an alternative boundary condition instead of (^. Under this new condi- 
tion we can repeat our analysis of the stationary solutions: Typically Eq = g{J) would now 
be different from a fixed point of the discrete mapping (pO]), so that the uniform field profile 
would not be a stationary solution of the discrete equations. The closest to the uniform 
solution would be a monotone field profile where a fixed point of ( pO]) is reached after a few 
iterations. For very long superlattices the difference with the results obtained with (|^) is 
small. Non-uniform two-domain solutions would still be possible: after a few iterates Ej 
would be close to the first fixed point of the mapping and the rest of our analysis in Section 
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H can be used to construct the non-uniform solutions. The differences would of course be 
more noticeable for short superlattices. If (for example) Eq = g{J) turns out to be small 
(think for instance of a linear law g{J) = rJ with a small resistivity r), there would be a 
region of the superlattice with a field smaller than that corresponding to the low field do- 
main of Section y. These differences should be experimentally testable, and in the absence 
of such tests and given the present good agreement with known facts , we prefer to stick 
to the simpler situation provided by the boundary condition (|]). Would future data force 
another choice of boundary condition upon us, our methodology could obviously be used to 
analyze the corresponding problem. 

Our theory compares favorably with previous ones: 

Grahn et al considered a continuous Poisson equation that ignored the holes and a 
stationary Ampere's law with a drift velocity consisting of a sum of delta functions. They 
imposed a step-like stationary electric field profile as a solution and obtained several qualita- 
tive features of the experimentally observed domains. Their results share with ours the fact 
that the average of the stationary current density at the branches representing non-uniform 
states is flat (independent of the applied bias. Section II). Dynamics and stability properties 
of non-uniform profiles were not considered in Ref. 0. 

Laikhtman's theory |jl2| consists of a discrete Ampere's law plus the condition that the 



total voltage bias is constant. He did not consider the carrier densities and equations for them 
(Poisson's law and a rate equation for the holes). Thus the only coupling between quantum 
wells in his theory comes from the bias condition and the electric field in the wells may 
alternate between values at the different branches of v{E) without restrictions. Then quite 
complicated non-monotone stable stationary field profiles with several domains are possible. 



as indicated in Ref. |T2[ . Within Laikhtman's theory these profiles cannot be eliminated: the 
only possible correlation between domains would be caused by assuming coherent tunneling 
transport between wells. Our model yields monotone non-uniform stationary profiles in a 
natural way without appealing to coherent tunneling. The number of possible non-uniform 
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stationary states is also much smaller than in Ref. . 

Korotkov et al's theory |]13| corresponds to quite different superlattices (slim) where 
the tunneling processes are single-electron tunneling events and the charge is quantized. 
Interesting effects are predicted in this limit: coexistence of Bloch oscillations and oscillations 
due to single-electron tunneling 



Other theories are based upon drift-diffusion models and are unable to reproduce im- 
portant qualitative features present in the experiments such as the multiplicity of steady 
states in the NDR region, and the eventual damping of the oscillations with time In 
fact, our model resembles known drift-diffusion models in the continuum limit (such as the 
one in Ref. |2^) for which the uniqueness of the steady state can be proved easily, at least 
in the limit of small diffusion coefficient (these continuum drift- diffusion models are also 



known to have Gunn effect oscillations p7| , p8| among their possible solutions |2^, which is 
not the case for Eqs. (p!2|-|16D). Thus it seems that the discreteness of our model is a crucial 
ingredient in reproducing both important static and dynamical features of transport in a 
superlattice made out of weakly coupled quantum wells. Work on derivation of discrete drift 
models from microscopic ones is now in progress. 
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APPENDIX: LINEAR STABILITY OF THE STEADY SOLUTIONS 



Linearization of Eqs. (P"^|-[T^) about any steady solution (denoted by Ej, n'j, p'j and J) 
yields the following system for the disturbances {ij, hj, pj and J) from it: 

Cj - = Uj - Pj , (Al) 
(3^ = J-v{E^)n,-n'^v'{E^)e,, (A2) 
^ = -n^P,-p%, (A3) 

Ee,=0, (A4) 
Co = ei . (A5) 

where, as usual, j runs from 1 to N. The system above can be reduced: we can eliminate the 
n's in favor of the p's and the e's by using (|A1|) and (|^). Moreover, summing ( [X^ ) for all 
j = 1, . . . , A^, and using (X|), we obtain J in terms of the same variables. Finally a system of 



2N equations (|X^- |A3| ) for the 2N unknowns {ij, pj, j = 1, . . . , N) results. Time evolution 
preserves a constant value for the sum of all the fields, but the condition that this value 
be zero is to be imposed as an additional constraint to the resulting system of equations. 
They can be written in matrix form as dx/dt = Mx, where we have defined the vector 
X = (ei, . . . , eN,pi, ■ ■ ■ ,Pn)^ ( ^ denotes the transpose). The eigenvalues of M should in 
principle determine the linear stability of the steady solutions; in practice, however, one of 
the eigenvalues of M is spurious and we should eliminate it by taking into consideration the 
constraint (|X^). 

We can find the exact eigenvalues for the uniform solution (p!^-[T7|). In this case, the 
matrix M is 



M 



(a a ^ 

^ee ^ep 



A A 



(A6) 



where 
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(3N p ' 



App = -2^7/ . 



(A9) 
(AlO) 



Here v and v' denote v{(f)) and v'{(f)) respectively, / is the N x N identity matrix and A, B 
and F are the foUowing N x N matrices: 

'^1 ••• -1^ 



A^ 



yl ••• -ly 



(All) 



B 





-11 ••• 

-1 1 ■-. : 
; ■-. ■-. ■-. 

••• -11 



V 



1 ••• 1 

The matrix M has only four distinct eigenvalues, given by the expressions: 

Ao = 0, 

Ai = -2V7, 



± 



(A12) 



(A13) 



(A14) 
(A15) 



(A16) 
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The two first eigenvalues are simple: Aq is the spurious eigenvalue, corresponding to the 
eigenvector (1, . . . , 1, 0, . . . , 0)"'" (notice that J2jLi % 7^ for this eigenvector) and therefore 
it has to be ignored; Ai corresponds to the eigenvector (0, . . . , 0, 1, . . . , 1)'^ and it is always 
real and negative. The remaining two eigenvalues (A+ and A_) have algebraic multiplicity 
— 1 each. It is not difficult to check that they are always real and that at least one is 
always negative. Accordingly, the stability of the uniform solution depends on the sign of 
the other eigenvalue. It is simpler to consider the sign of the product A+A_. The solution 
will be unstable whenever this product is negative and neutrally stable when it is zero, i.e. 
when 

V + 2^v' = . 

This equation is exactly p^)! Its minimum is 71, the photoexcitation below which the 
uniform solution is the only existing stationary state. 

It is clear that we could in principle ascertain the linear stability of any steady solution of 



(|IJ-|T6|) by the method so far explained. We would use any of the well known diagonalization 
numerical routines on the matrix M corresponding to the steady solution. However these 
routines quite disastrously fail to determine the exact value of the eigenvalue of M because 
of the high multiplicity thereof: the errors of the numerically determined eigenvalues can 
be as high as 10%, and, what is even worse, they become artificially simple and sometimes 
even acquire a spurious imaginary part. The fact that the matrix M for the uniform steady 
solution is not diagonalizable increases the numerical difficulty. It is important to notice 
that the linear stability of the uniform stationary solution is not affected by these numerical 
errors: it depends only on the sign of the largest eigenvalue and we have checked that this 
sign is given correctly by the numerical code (we have always compared the results provided 
by the code with the formula ( [A16| ) with the + sign that yields the largest eigenvalue 
exactly). 

It could seem that these numerical troubles due to eigenvalue degeneracy may not appear 
for domain-type steady solutions (which are non- uniform) . However, we do not have now 
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exact means of checking this conjecture at our disposal Thus our numerical determination 
of the hnear stabihty of domain-type steady solutions has always been confronted with the 
results of direct simulation of the full system of equations whenever the comparison was 
possible. We have found that (as in the case of the uniform solution) the linear stability 
of the non-uniform steady states is correctly predicted by the numerical determination of 
the eigenvalues. There is a reasonable agreement between the two numerical procedures for 
the real part of the largest eigenvalue while there are large errors for the imaginary part. 
Thus while the eigenvalue calculation gives acceptable estimations of the damping of the 
photocurrent oscillations, it provides poor estimations of the frequency. Hence the stability 
properties of a given stationary state are well determined by the numerical procedure and 
we have used it, together with the simulations, to delimit the region of the 7 — phase dia- 
gram (Fig. 1^) inside which stable non- uniform stationary solutions exist. Work on devising 
reliable numerical codes to determine the oscillation frequency from the eigenvalues is now 
in progress. 
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FIGURES 

FIG. 1. Different regions and processes that contribute to the PL spectrum (see text). Only 
the conduction band is represented. 

FIG. 2. Velocity of the electrons as a function of the electric field. It is modeled by two 
lorentzians centered at Em and E^' which correspond to the sequential resonant tunneling current 
between the quantum wells j and j + 1. At a field Em, there is SRT between the first level of the 
j-th well and the second level of the (j + l)-th well, whereas at a field Em' there is SRT between 
the second level of the j-th well and the third level of the {j + l)-th well. (In different samples, the 
latter SRT process may relate the first level of the j-th well and the third level of the (j + 1)— th 
weh). 

FIG. 3. Static characteristic curve: current ( J) versus applied bias (</») for laser power 7 = 0.016 
(/?, which is only relevant for the stability of the different solutions - see ( [A2| ) in the Appendix, is 
assumed to be 0.05). The solid line corresponds to the stable uniform solution. In the NDR part of 
the curve there is a region where this solution becomes unstable (short-dashed line) . The remaining 
curves correspond to two-domain solutions. There are 39 branches (one less than quantum wells), 
though only one every three branches has been plotted for the sake of clarity. The inset enlarges a 
small region of the figure showing the branches in full detail. Each branch has a different number 
of wells (from 1 to 39) in the low field domain (the number written on top of some of them) . Every 



branch consists of two curves, according to the way they are built with the discrete mapping (20) 
(see text for details): one plotted with triangles, corresponding to jumps from the first to the third 
fixed points of ([20|), and another plotted with a combination of squares (stable part) and long 
dashed lines (unstable part), corresponding to jumps from the second to the third fixed points of 
(pOl). Important features to stress from this figure are: first, the coexistence, for a given (p, not 
only of several domain solutions but also of domain solutions and the uniform solution (see also 
Fig. ^), something relevant to the appearance of hysteresis and memory effects, and second, the 
average "flatness" of J in the bias region between the two local maxima of v{E). 
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FIG. 4. The discrete mapping: f{E; 7, J) and E versus £^ for 7 = 0.016 and five different values 
of J. (a) For small J (below the minimum of y/jv{E)), only one fixed point of the discrete mapping 
appears; thus only uniform solutions are allowed, (b) There are three fixed points of the discrete 
mapping, but only 2 — 3 jumps are allowed, (c) Now both 1^3 and 2^3 jumps are allowed. 
The continuous lines show the discrete mapping process used to construct a {1 ^ 3) non- uniform 
stationary solution with two domains: After remaining at the low field domain value E^, the field 
jumps at the (m+ l)-th well to a different solution of f{E;j, J) = El. In successive iterations the 
electric field tends to its high-field value Eh- The short-dashed lines show the same construction 
for a 2 — > 3 solution, whereas the long-dashed lines correspond to a different (unstable) type of 
2 — 3 solution explained in the text, (d) Same as (c), but showing the coalescence of solutions 
with 1 — 3 and 2 — 3 jumps for J = y/jv{EM)- (e) For large J, above the first positive maximum 
of ^v{E), there is only one fixed point of the discrete mapping and the situation is the same as 
in (a): the only possible stationary solutions are uniform. 

FIG. 5. Phase diagram of the model: Laser power (7) vs. voltage bias {<p). Stable non- uniform 
solutions with domains exist only inside the continuous curve. Inside the dashed line the uniform 

solution is unstable. The shaded rectangle approximately delimits the region where oscillations of 
the PC (i.e., of the domain wall of the electric field profile) can be obtained. 

FIG. 6. Evolution of the domain wall in time for 7 = 2 x 10~^ and ^ = 1.2. (a) is a contour 
plot with time flowing in the y-ax;is and the quantum well number in the x-axis. The electric 
field contours go from 0.5 to 1.7 in 0.1 steps. At t = wc start with the uniform solution. The 
perturbation is clearly visible at the first well at short times. The creation and the displacement of 
the domain wall can be visualized. At t = 0.5 the two domains are already defined (the darker the 
figure the higher the electric field) and the domain wall begins to oscillate. Most of the wells have 
an electric field centered around the first maximum of v{E) (low field domain). In (b) we represent 
the electric field profile at short times (domain formation) and in (c) the domain wall oscillations. 
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FIG. 7. The PC (J) evolution for the same values of the parameters as Fig. ^ The transient 
before the oscillations corresponds to the formation of the low field domain. The oscillations of the 
domain wall induce those of the PC. The maxima are reached when the domain wall is sharper. 
In the inset it can be observed the small damping of the PC, as the system slowly approaches the 
stable uniform solution. This damping disappears for large enough laser power (see text). 

FIG. 8. Evolution of the PC, J, for different values of the photogeneration, 7 (in units of 10^^), 
which increases from bottom to top of the figure, as indicated at the right margin. The voltage 
{(j) = 1.1) is the same for all the curves. 

FIG. 9. (a) Electric field profile and (b) charge at each well corresponding to a maximum (full 
line) and a minimum (dashed line) of the PC. 

FIG. 10. Oscillations of the PL spectrum corresponding to Fig. |^ 



34 



